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Abstract The superconducting and structural properties of S/F/S (Super- 
conductor/Fcrromagnct/Supcrconductor) heterostructures have been studied 
by means of microwave measurements (1-20 GHz) and x-ray absorption fine 
structure (XAFS) spectroscopy. Nb/PdNi/Nb trilayers have been studied as 
a function of F layer thicness. With respect to pure Nb, XAFS analysis shows 
that the heterostructures exhibit larger structural disorder in the S layers. Mi- 
crowave measurements show evidence for a progressively weaker vortex pinning 
with increasing F thickness. However, no clear correlation is found with the 
local disorder in Nb: the weakest pinning is not in the most disordered trilayer. 
Therefore the structural disorder in the superconducting material cannot ex- 
plain on its own the changes in vortex pinning. We argue that the F layer acts 
on the superconducting state itself. We propose possible explanations for the 
observed behavior. 
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1 Introduction 

The physics of Superconducting/Ferromagnetic (S/F) heterostructures allows 
to investigate the competing magnetic (F) and superconducting (S) ordering 
in artificial structures where a spatial separation among S and F exists. The 
interest in such systems has been revamped in recent years [2[2] , in particular 
with respect to the many effects that can arise when the ferromagnetic coher- 
ence length is of the order of the F layer thickness. Similarly, a new interest in 
the vortex physics (namely, vortex pinning) in S/F structures has arosen pH 

®. 

However, measurements of the vortex state response at high frequencies 
(radio and microwave) are lacking, despite the important information that 
they can give, as it was recognized very early [5]. In particular, high frequency 
measurements can give direct measures of several vortex parameters, such 
as the pinning constant (or the depinning frequency) and the free-flux-flow 
resistivity (see below, Sec|3]). 

In this paper, we present a preliminary investigation of the pinning prop- 
erties of S/F/S trilayers (and pure Nb as a comparison) using swept-frequency 
microwave measurements in the frequency range 1-20 GHz. Since defects in 
the material can affect severely the vortex pinning, we also performed mea- 
surements of the local structure of Nb in our Nb/PdNi/Nb trilayers in order 
to unveil possible correlations between pinning properties and local disorder. 
To this aim we exploited the Nb-K edge XAFS technique. 

The paper is organized as follows: in Section [2] we shortly describe sample 
preparation and we report measurements of the local structure around Nb 
by means of XAFS. In Section [3] we report and discuss our measurements of 
the vortex parameters obtained from swept-frequency microwave response. In 
Section [4] we summarize the results. 



2 Samples and local structure measurements. 

Pure Nb films and Nb/Pds4Nii6/Nb S/F/S trilayers have been grown by ultra- 
high-vacuum, dc magnetron sputtering. All samples were grown at room tem- 
perature. Extensive details on sample preparations have been given elsewhere 
[5j. Single-crystal sapphire substrates were used, in order to avoid the compli- 
cations arising with the analysis of the microwave response when Si substrates 
are employed [91IT0]. Three heterostructures were studied. They consisted of 
two Nb layers of nominal thickness ds= 15 nm, and an inner PdNi layer with 
thicknesses dp=l, 2 and 8 nm. A pure Nb sample of nominal thickness djvb=30 
nm has been examined for comparison. The thickness of the samples has been 
carefully calibrated against the deposition rate. As expected, the critical tem- 
perature T c dropped with increasing d,F [H], and we obtained T c = 7.5, 6.2, 
5.1, 4.1 K in the samples with dp= (pure Nb), 1, 2, 8 nm, respectively. 

In this paper we are mostly interested in the possible effects of the fer- 
romagnetic thickness on the vortex state properties, which in turn might be 



Superconducting and structural properties of Nb/PdNi/Nb trilayers. 



3 



sensitive to the structural disorder. Local disorder can affect the supercon- 
ducting properties in several ways: it can act as a pinning center, thus giving 
rise to e.g. large critical currents; it can weaken the stiffness of the superfluid, 
thus giving rise to a smaller superfluid density; it can increase the quasiparti- 
cle scattering rate, thus affecting the normal state resistivity and possibly the 
flux-flow resistivity. Thus, an evaluation of the local disorder is certainly an 
essential information. To this aim we exploited the XAFS (x-ray absorption 
fine structure) technique to probe the local atomic structure in the pure Nb 
film and in heterostructures with dp=2 and 8 nm . 

Nb-K edge (£jvf,=18.986 keV) x-ray absorption spectra (XAS) were col- 
lected at the GILDA-BM08 beamline [T^] at the European Synchrotron Ra- 
diation Facility (ESRF Grenoble, France) in fluorescence geometry. GILDA 
beamline optics employs double independent crystals Si [3 11] monochromator 
and a couple of Pd-coated mirrors for harmonic rejection. Focused x-ray spot 
(approximately 0.2mm (v) x 1 mm(h)) is obtained by vertically bent mirrors 
and sagittally focused second monochromator crystal. The x-ray absorption 
spectra were collected in fluorescence geometry keeping the samples at 77 K 
to reduce the structural disorder induced by thermal vibrations. The sam- 
ple surfaces were oriented about 45° from the polarization direction of the 
incoming x-ray beam in order to reduce the polarization effects in the spectra. 

Nb k a fluorescence intensity [Ek„ ~ 16.6 keV) was measured using a 
13-element ultrapure Ge multidetector. Bragg reflections were identified by 
comparing the spectra collected with the different elements and removed 
by interpolating the data under the peaks; the total fluorescence signal is cal- 
culated as If = The incident intensity Iq was measured using an Ar 
filled ionization chamber to calculate a(E) = If/Io, proportional to the Nb 
linear absorption coefficient fi(E) (the thin Nb films ensure negligible self ab- 
sorption effects). Up to 8 spectra were measured for each sample and averaged 
up to improve the data statistics. The average statistical noise on the data, 
calculated as described in (T3], results in the range 1.5 — 1.8 • 10~ 3 (Table 1). 

The ESTRA and FitEXA codes [T3] were used for XAS treatment and 
XAFS quantitative analysis: the a(E) signals were treated using standard pro- 
cedures for background subtraction, normalization and extraction of x-ray ab- 
sorption fine structure (XAFS) signal x exp (k) [13] (k = 7T 1 yj2m e {E - E Nb ) 
is the photoelectron wavenumber). 

Quantitative XAFS data analysis has been performed by least square data 
refinement procedure in the reciprocal space, by fitting the raw A:-weighted 
spectra k\ exp to the theoretical XAFS formula [M] : 

k x th {k) = S 2 Q J2 Bin(2Jfcrj + <Me"W^ (1) 

Eq. ([I]) is a sum of partial contributions coming from single and multiple scat- 
tering terms which are calculated assuming a Gaussian model for the pair (or 
multiple scattering) distribution functions. Nj, rj and cr| are the multiplicity 
(coordination) numbers, interatomic distances and mean square relative dis- 
placement (MSRD) of the atomic pair distribution. The S% is an empirical 
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r, = N T = 8 

rn = a N n = 

r rn =as/2 N HI =12 
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Fig. 1 Local geometry in bee structure: the five coordination shells used in the analysis 
are highlighted. The bec geometry constrains the scattering path lengths r{ and multiplicity 
numbers 2Vj of the coordination shells used for XAFS data refinement as summarized in the 
table. 



parameter accounting for many body losses, and it was fixed to 0.8 for all the 
samples. The A j and 4> j (amplitude and scattering functions) and X j (the 
photoelectron mean free path) were calculated using the FEFF 8.2 code [15] . 
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Fig. 2 Experimental Nb K edge XAFS data (points) and best fit (full line) for all the 
investigated samples (vertically shifted for the sake of clarity). 



Nb has body centered cubic (bec) crystallographic structure, depicted in 
Figjl] with lattice parameter (cube edge) a=3.30A. The five coordination shells 
shown in Fig. 1 were used to refine the Nb K edge XAFS data. Notice that 
the aligned Nbo— Nb/— Nby configurations along the bec cube diagonal give 
rise to strong multiple scattering contributions which have been considered 
in the analysis. The multiple shell data refinement has demonstrated to im- 
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Sample 


Ri 


I 


<4 


AE 


< 


R 2 




[A] 


[x A 2 ] 


[x A 2 ] 


[eV] 


[xicr 3 ] 




Nb(30) 


2.85(1) 


5.4(4) 


5.8(6) 


-3.9 


1.8 


0.13 


Nb(15)PdNi(2)Nb(15) 


2.84(1) 


7.7(6) 


8.7(7) 


-4.9 


1.7 


0.11 


Nb(15)PdNi(8)Nb(15) 


2.85(1) 


7.3(6) 


6.9(7) 


-4.2 


1.6 


0.12 



Table 1 Best fit results of Nb K edge XAFS data analysis. 



prove the XAFS data analysis reliability [15111 7| . as it allows to apply severe 
constraints among the structural parameters. In order to reduce the number 
of free parameters and the correlations among them all the 77 distances were 
constrained to the lattice parameter a (Fig.l) and the multiplicity numbers 
were fixed to their crystallographic values. The MSRD and the energy scale 
shift (AE) were refined. An example of best fit, showing the partial contribu- 
tions and the residual, is reported in Fig(3J The experimental data and best 
fit for all the three samples investigated are also shown in Figj2] The squared 
residual function: i? 2 [13] is reported (Table 1) to evaluate the best fit quality. 




j » s 5 (SS+MS) 

Residual - 
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Fig. 3 Example of XAFS data analysis, Nb(30) sample: experimental k\ exp {k) XAFS data 
(points) and best fit (full line) are shown at the top. The partial contributions (coordination 
shells) are shown below, vertically shifted for sake of clarity, labeled according to the bee 
geometry (FigJlJ. The residual (experimental minus best fit) is shown at the bottom. 

The best fit parameters provided quantitative details about the Nb lo- 
cal structure (Table 1): firstly we notice that the geometrical constraints 
based on bec structure provide a satisfactory best fit in Nb(30) as well as in 
Nb(15)/PdNi(x)/Nb(15) trilayers, thus confirming the validity of the model. 
The lattice parameter corresponds to what expected in bulk Nb. Therefore, the 
Nb local structure is close to the bulk Nb structure in all samples, and only 
weakly affected by the magnetic PdNi layer. Nevertheless, the MSRD (and 
thus the structural disorder) is larger in trilayers (about 40% higher) than in 
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pure Nb. The growth of the interlayer has then some effect on disordering the 
Nb atomic structure at the local scale. Interestingly, this effect is nearly in- 
dependent from the interlayer thickness, as Nb(15)/PdNi(8)/Nb(15) exhibits 
comparable, or even smaller, disorder with respect to Nb(15)/PdNi(2)/Nb(15). 

3 Microwave measurements. 

We measured the real part of the complex microwave resistivity in the wide 
frequency range 1-20 GHz, by means of the so-called Corbino disk technique. 
The method has been extensively described elsewhere [18] , In short, the thin 
film terminates a coaxial line. From measurements of the complex reflection 
coefficient of the line it is possibile, after suitable calibration, to obtain the 
frequency dependence of the real part of the resistivity. Measurements are 
performed at fixed field H and temperature T, and the frequency is swept 
up to 30 GHz. The noise increases significantly above 20 GHz, so that in our 
calibration setup the measurements are reliable only up to ~ 20 GHz. We 
present here a set of homogeneous measurements, taken in samples with dp 
= 1, 2, 8 nm, and in the pure Nb film, at the same reduced temperature 
t = T/T c ~0.86 and at the same reduced field h = H/H c2 ~0.5. 

The frequency dependent resistivity allows for an accurate determination 
of the vortex parameters. On one side, the microwave field induces very short- 
range oscillations in the vortex system. The mean-square displacement of a 
vortex line from its pinning site depends on the microwave frequency and field 
strength, but it is typically much less than 1 nm |20j . This peculiar feature 
makes the analysis of the data quite straightforward. In particular, the vortex 
dynamic can be assumed to be "single- vortex-like" : since the displacement 
of each vortex is much smaller than the intervortex separation, which with 
the flux density here applied (between 0.125 T and 0.25 T) is in the ^100 
nm range, no perturbation of the static configuration is determined by the 
microwave current. In this case, the generalized model for vortex motion |21j . 
which recovers many well-known models [5j l22|l23] in the appropriate limits, 
can be applied. The vortex motion complex resistivity reads: 

Pvm — Pvm.l + l Pvm.2 — P f f ~. — ^ 1 \ Z ) 

l + iv/vo 

where the flux-flow resistivity pff has been discussed previously [TS], is 
a characteristic frequency separating the pinned regime (y < vq) from the 
flux-flow regime (v > vq), and the creep factor e < 1 takes into account 
thermally activated processes. The so-called depinning frequency v p , which is 
a commonly used parameter to quantify the strength of vortex pinning, can 
be easily derived from v$ and e, once a specific model is assumed to obtain 
the analytical expression for e. All the details have been discussed thoroughly 
elsewhere [21] . 

It is worth mentioning that all three vortex parameters can be measured in- 
dependently even from the real part of the frequency-dependent vortex motion 
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resistivity. Typically, curves of p V m,i(y) are fitted to Eq.([2|, yielding j/o, e. 
Here, in order to derive v p , we adopted the Coffey-Clem model [22], but we 
emphasize that different models, such as [53], gave basically the same results 
within a scale factor. 

In most cases [24j the microwave complex resistivity has been obtained at a 
single frequency, typically with the help of some kind of resonating technique, 
and the analysis of the data has relied on the assumption of vanishing flux 
creep, which we show below is not the case in our measurements. Thus, in the 
present case, the frequency-dependent measurements are invaluable in that 
they allow for the simultaneous measurements of all three vortex parameters. 

In Figj4] we report typical swept-frequency data taken in one of our sam- 
ples, together with the fit obtained from Eq.([2|. Apart from the oscillations, 
which have been described in [25], the overall fit is quite satisfactory and 
yields the vortex parameters with a small uncertainty. The vortex parame- 




5 10 15 v(GHz) 

Fig. 4 Full dots: normalized vortex motion resistivity, p vm ,\/pn vs frequency; continuous 
line: fit with Eq.Q. To avoid crowding, only 15% of the data are shown. Oscillations are 
typical of sudden increase of dissipation, that has been discussed elsewhere |25| . Sample is 
Nb/PdNi(2nm)/Nb. T=4.4 K, HoH= 0.2 T. 



ters are reported in Figj5j here we observe that the ferromagnetic thickness 
cLf has a dramatic influence on pinning and creep. In pure Nb the depinning 
frequency v v «3.5 GHz points to rather strong pinning, in agreement with 
other estimates [26] , The sample with thinner dp does not change much in 
pinning, but the sample with dp=2 nm shows an abrubt collapse of the de- 
pinning frequency to v v w 2 GHz. By contrast, thermal activation increases 
sharply immediately when an F layer exists, and steadily keeps increasing with 
increasing dp. Thus, S/F/S multilayers clearly show a much weaker pinning 
and larger thermal activation. 

By comparing the results of the present Section with the local structure 
measurements of Secj2] we come to the interesting conclusion that if a corre- 
lation between the disorder and the change in v p and e exists, it is far from 
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being trivial. In fact, while the local structural disorder is larger in S/F/S 
trilayers than in the pure Nb film, the trend is not a monotonous function of 
dp. The larger disorder is found by XAFS in the sample at dp— 2 nm, which 
does not exhibit neither the weakest pinning, nor the larger thermal activa- 
tion. Thus, we must argue that the magnetic properties of the F layer play a 
role in modifying the superconducting response of trilayers. 

In fact, both the depinning frequency and the creep factor are related to 
the condensation energy: the larger the condensation energy of the condensate, 
the larger the pinning strength, the smaller the creep factor. Thus, it can be 
reasonably argued that the F layer affects vortex pinning and creep through a 
reduction of the condensation energy. This explanation would be in agreement 
with the observed reduction of the superfluid fraction observed in S/F bilayers 




2 d F ( nlTI ) 6 8 

Fig. 5 Pinning parameter u p (depinning frequency), right axis, and creep parameter e, left 
axis, as a function of the ferromagnetic layer thickness dp. As it can be seen, the increase 
of the ferromagnetic thickness determines a decrease of pinning strength and an increase of 
thermal creep. All measurements have been taken at T/T c ~0.86 and H/H c 2 ~0.5. Lines 
are guide to the eye. 



4 Summary. 

We have presented measurements of the frequency-dependent vortex motion 
resistivity in pure Nb and Nb/PdNi/Nb trilayers. We have shown that vortex 
pinning and vortex creep arc heavily affected by the presence and thickness of 
the F layer: pinning decreases and creep increases with increasing F thickness. 
Details about the local atomic structure in S phase (Nb) obtained from XAFS 
analysis demonstrate that structural disorder has an effect in modifying the 
pinning and creep strength, however the structural disorder alone cannot ex- 
plain the observed weakening of pinning/strengthening of creep as a function 
of F-layer thickness. Combined XAFS and microwave results are relevant to 
recognize the role of F layer in modifying the superconducting state itself. 
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While deeper investigations are necessary to clarify fully the effect of the F 
layer on the superconducting properties of the heterostructures, we argue that 
the effects observed in the vortex response might arise from a reduction of the 
condensation energy. 
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